Thermodynamic signatures for topological phase transitions to Majorana and Weyl superfluids in 
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We discuss the thermodynamic signatures for the topological phase transitions into Majorana and Weyl super- 
fluid phases in ultracold Fermi gases in two and three dimensions in the presence of Rashba spin-orbit coupling 
and a Zeeman field. We analyze the thermodynamic properties exhibiting the distinct nature of the topological 
phase transitions linked with the Majorana fermions (2D Fermi gas) and Weyl fermions (3D Fermi gas) which 
can be observed experimentally, including pressure, chemical potential, isothermal compressibility, entropy, 
and specific heat, as a function of the interaction and the Zeeman field at both zero and finite temperatures. We 
conclude that among the various thermodynamic quantities, the isothermal compressibility and the chemical 
potential as a function of the artificial Zeeman field have the strongest signatures of the topological transitions 
in both two and three dimensions. 
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I. INTRODUCTION 



Ultracold fermionic quantum gases have many advantages 
for investigations of quantum phenomena in a highly tunable 
way. The experimental ability to control the interaction be- 
tween the Fermi atoms via Feshbach resonances [ 1 ] enables to 
study the crossover phenomena [j2|] from the BCS superfluid 
of weakly bound Cooper pairs to the BEC of tightly bound 
molecules |3], universal properties in the unitary regime [Ql, 
vortices |5], and so forth. In addition, the experimental con- 
trol of the internal spin states via radio-frequencies provides 
the study of the emergence of new superfluid phases in a 
Fermi gas with population imbalance which serves as an ef- 
fective Zeeman field [6- 8]. Another new development enables 
to trap two dimensional fermions in the quantum degenerate 
regime lUIH. 

Furthermore, in recent experimental breakthroughs the 
spin-orbit coupling (SOC) has been developed using two- 
photon Raman processes in bosonic as well as fermionic 
gases lITlUrsll . which has given rise to theoretical investiga- 
tions on the new quantum states of matter lfT6l - [l9ll . Thus far, 
the experimental realization of SOC is an equal combination 
of Rashba li2Ql and DresseUiaus [21] types of SOC. It has 
been shown theoretically that in the presence of Rashba SOC 
and a perpendicular Zeeman field, 2D cold fermion systems 
with an i-wave Feshbach resonance can enter into a topo- 
logical superfluid phase with zero energy Majorana fermion 
excitations localized at the order parameter defects such as 
vortices ifTsl l22l l23ll . It has also been shown that in 3D the 
corresponding system is capable of supporting the intriguing 
Weyl fermion excitations [17]. The realization of ultra-cold 
fermionic gases with Rashba SOC and a Zeeman field thus of- 
fers the tantalizing possibility of experimental engineering of 
all three kinds of fundamental relativistic fermions, namely, 
the Dirac, Weyl, and Majorana fermions llbil . In addition 
to these powerful experimental and theoretical developments, 
there have been more recently added to the list. One is the 
development of the methods for high precision measurements 
of the thermodynamic properties such as the compressibility. 



the chemical potential, the entropy, and the heat capacity for 
uniform ^-wave Fermi superfluids |25, 26]. In recent devel- 
opments, Ref. li26ll has completely determined the universal 
thermodynamics of the interacting Fermi gas across the s- 
wave superfluid transition and has captured the onset of su- 
perlfuidity in the compressibility, chemical potential, entropy 
and the heat capacity. The question we address in this paper is 
what are the corresponding thermodynamic signatures for in- 
teracting ultra-cold fermions across the topological superfluid 
transitions to phases supporting Majorana fermions (2D) and 
Weyl fermions (3D). Our main conclusion is that among the 
various thermodynamic quantities, the isothermal compress- 
ibility and the chemical potential as a function of the artifi- 
cial Zeeman field have the strongest signatures of the topo- 
logical transitions. The signatures for the transition from the 
non-topological i-wave superfluid to the Majorana superfluid 
phase in 2D are stronger than the corresponding signatures in 
3D where the transition is in the Weyl superfluid phase. The 
other thermodynamic quantities such as the pressure, entropy, 
and specific heat have only weak signatures of the topological 
transitions (that is, the transition between a non-topological 
superfluid and a topological superfluid as a function of the 
Zeeman field), even though they are useful to capture the topo- 
logical superfluid to the normal phase (pair potential A = 0) 
transitions in both two and three dimensions. 

In condensed matter systems, SOC, which results from 
the coupling between the spin degrees of freedom and the 
orbital motion of a quantum particle, plays an important 
role in understanding interesting physics associated with the 
topological quantum phenomena including the spin Hall ef- 
fect 11271 l28ll . topological insulators |29], and topological su- 
perconductors ||30ll- These phases generically arise as a result 
of topological quantum phase transitions (TQPT) in the ap- 
propriate systems. TQPTs separate distinct phases of matter 
which have exactly the same symmetries and thus are not as- 
sociated with any spontaneous symmetry breaking in the un- 
derlying Hamiltonian. Yet, the phases separated by TQPTs 
have distinct topological properties which can often manifest 
themselves in the properties of the order parameter defects . 
For instance, it has been pointed out recently that spin-orbit 
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coupled electron- or hole-doped semiconductor thin films 
(2D) or nanowires (ID) with proximity induced i-wave su- 
perconductivity and a suitably directed Zeeman field (F) sup- 
ports novel non-Abelian topological states with the order pa- 
rameter defects supporting localized to pologic al zero-energy 
excitations called Majorana fermions ll3ll - l36ll above a criti- 
cal Zeeman field (F > F^). These theoretical developments 
in the condensed matter system followed earlier similar pro- 
posals in the context of ultra-cold fermions in two dimensions 
in the presence of artificial SOC and a Zeeman field and with 
i-wave pairing interactions ll22ll . While the two dimensional 
cold fermion systems support the Majorana fermions in the or- 
der parameter defects (but are otherwise completely gapped in 
the bulk), the corresponding three dimensional systems have 
topologically protected gapless nodes in the bulk supporting 
the relativistic Weyl fermion excitations. With Rashba-type 
of SOC and an increasing strength of the Zeeman field, the 
three dimensional system supports a series of 3D topological 
quantum phase transitions |37] from a non-topological super- 
fluid state with fully gapped fermionic excitations to a topo- 
logical superfluid state with four protected Fermi points (i.e., 
gap nodes) and then to a second topological superfluid state 
with only two protected Fermi points 1 17|. Such 3D topologi- 
cal gapless superfluids with Weyl fermions, which result from 
quasiparticle excitations in a linearized Hamiltonian near the 
nodes, is different from the 2D fully gapped topological su- 
perfluids with Majorana fermions [18,.31J. Recently, similar 
Weyl fermion physics in an analogous 3D condensed matter 
system has also been discussed in the context of ferromag- 
netic superconductors |38]. The results on the thermody- 
namic signatures of the topological transitions to Majorana 
and Weyl superfluid phases discussed in this paper apply to 
both cold fermion and condensed matter systems although our 
focus will mostly be ultra-cold fermions where the methods 
for high precision measurements of the thermodynamic quan- 
tities have recently been successfully developed. 

The paper is organized as follows. In Sec. II, we derive 
the mean-field order parameter and atom-number equations 
for ultra-cold fermions with Rashba SOC and a Zeeman field. 
We present the zero temperature phase diagram and the ther- 
modynamic quantities as a function of interaction for both 2D 
and 3D Fermi superfluids from the self-consistent numerical 
calculations in Sec. III. Section IV provides the phase dia- 
gram and the corresponding thermodynamic properties at fi- 
nite temperature. Section V consists of discussion and con- 
clusions. 



the single particle Hamiltonian 



yy' 



(1) 

where y =t, X are the pseudo-spin of the atoms, ^ = 
-fi^V^/2m - yu, /i is the chemical potential, a is the Rashba 
SOC strength, / is the 2x2 unit matrix, cr, is the Pauli matrix, 
F is the strength of the Zeeman field, and c^y is the atom an- 
nihilation operator Him - -9 J dT^cl-^cl^CriCrt is the ^-wave 
scattering interaction with g = ATrH^aJm, and the scattering 
length as can be tuned by the Feshbach resonance. 

The grand partition function at temperature T is Z = 
J D[if/, if/] exp {-S [tf/, i^]) with action 
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Y^4,yir,T)—i(,y(r,T) + 'H(r,T) 



(2) 



Using the standard Hubbard-Stratanovich transformation that 
introduces the pamng potential A(r, r) - 0((Aj,(r, T)(/rf(r, r)) 
and integrating over the fermion variables lead to the effective 
action 
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where ^ = -H V /2m - fi and the matrix M is 
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dr+i^ -A 

-hi dr+li ^ 

A* dr-l] hi 

, -A* h^ dr-ll 
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where /i^ = -ifia{dy + idx) corresponds to the transverse com- 
ponent of the SOC field, F to the parallel component with re- 
spect to the quantization axix z, and f = ^ + F. 

To proceed, we use the saddle point approximation 
A(r,T) = Ao + ?7(r,T), and write M - Mq +Mp. The ma- 
trix Mo is obtained via the saddle point A(r, r) — > A which 
takes M Mo, and the functional matrix Mf = M - Mo de- 
pends only on rjir, t) and its Hermition conjugate. Thus, we 
write the effective action as5'etf = 5o+5/r. The first term is 
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where w„ 



(2n + \)nT . Here E^_^ are the eigenvalues of 



11. MEAN-FIELD THEORY 



The system we consider is a uniform i-wave fermionic su- 
perfluid with the atom density n - N/V in the presence of 
Rashba-type spin-orbit coupling (SOC) in the xy plane and a 
Zeeman field along the z direction. The dynamics of the Fermi 
gas can be described by the Hamiltonian H = Ho + Him, where 



Mk = 



f^k + r -Mk) -A 

hl(k) ^K-F A 

A* -fk-r /i^(-k) 

-A- hA-k) -^k + rj 



(6) 



which describes the Hamiltonian of elementary excitations in 
the 4-component Nambu spinor = (ck-f, Ckj,, c^,^^, c^|^p' . 
Here, /!j_(k) = a(ky + ik^) and the quasiparticle excitation en- 
ergy 



(7) 
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is the eigenvalue of Mk, where A = + correspond to the 
particle and hole branches, £0 = -^^r2(^TlAp)TcA^, 

and k± - + ky. For a - F - 0, Eq. (|7|i reduces to 

= A + |Ap in the standard BCS theoiy. 
The saddle -point grand potential function Q = -TlnZo is 



g 



(8) 



leading to the order parameter via minimization of Q with re- 
spect to A: 



A 1 , 



k,;;,/l 



dA* 



(9) 



where f(x) - 1/(1 + exp(-ji:/r)) is the Fermi-Dirac distribu- 
tion function. For the case of 2D gas, Eq. (|9]l can be expressed 
as 



2 {l-nr^/Eo)tmh(El^/2T)/4El^ 



1 



2ek + Eb 



where Ei, = h^/ma^ is the binding energy controlled by Fesh- 
bach resonances, and for 3D gas. 



Vm 
AntP-a 



= - Z [(1 - '/rV£o)tanh(£,-,^/2r)/4£,+„ - -L 



The ultra-violet divergence at the large k in Eq. (|9]) has been 
regularized through subtracting the term 1 /2ek and a, is de- 
fined as the renormalized scattering length. The total number 
of atoms can be obtained from N - - (dH/d/u): 



N-- 



^ Eq 



2E, 



tanh(£^/2r) 



(10) 

We self-consistently solve the order parameter equation 
(|9]l and the number equation (fTOl i for different parameters 
(aKf, r, V, T) for a fixed atom density n to determine A and jj. 

Here v - \ IKfGs and Kf - yin n\ is the Fermi momentum 
for a non-interacting 3D Fermi gas with the same density at 
Y - a - Q and v - Eb and Kf - {Inn)^!^ for a non-interacting 
2D Fermi gas. The energy unit is chosen as the Fermi energy 
Ep = n^Kj/2m. 



III. TOPOLOGICAL PHASE TRANSITIONS DURING 
BCS-BEC CROSSOVER (T = 0) 

A. Phase diagram 

In Fig. [T] we show the mean-field zero temperature phase 
diagram as a function of the Zeeman field F and the interac- 
tion parameter v in the presence of Rashba SOC {aKf - 0.4). 




1.5- 



(b) ' 


/ 






■ TS4 

— 1 1 


j NTS " 



0.5 1 1.5 2 
v=Eb/Ef 



-1 1 

v=1/KFas 



FIG. 1 . (Color online) Plot of (T - y) phase diagram for (a) 2D and 
(b) 3D Fermi gas at T = with aKf = 0.4. Here, v represents in- 
teraction parameters. Eg /Ef for 2D and 1 /KpOs for 3D systems. We 
label fully gapped non- topological superfluid phase as NTS, fully 
gapped topological superfluid phase as TSo, and topological super- 
fluids with 2 and 4 nodes as TS2 and TS4, respectively. Red and 
blue lines denote = -sjfi^ + and Fa = A, respectively, where the 
topological phase transitions occur. Green dashed lines correspond to 
fi = 0, below which /j > and above which fi < 0. Note that TS4 su- 
perfluid phase, which is characterized by 4 nodes along z-axis, arises 
when /J > and merges at the tricritical point, which is indicated by 
purple circle in the phase diagram of 3D fermion superfluid. 



We label the uniform topologically non-trivial i-wave super- 
fluid phases with zero, one, two, and four nodal points as TSo, 
TS 1 , TS2, and TS4, respectively. For 2D Fermi gas as in Fig. [T] 
(a), topologically non-trivial superfluid phase with one node, 
TSi, occurs only along F = Fc = -^Jfi^ + A^ line. Even though 
below and above F = Fc, the superfluid gap remains finite 
everywhere in the momentum space, thus the excitation spec- 
trum is fully gapped, F^ is a phase boundary between the non- 
topological superfluid (NTS) and non-trivial (TSo) superfluid 
phases. Even though the TSo phase is fully gapped, the phase 
is topologically non-trivial because in the presence of an order 
parameter defect such as a vortex, the vortex core supports a 
non-degenerate zero energy Majorana fermion excitation. The 
TSo phase has recently received wide attention in condensed 
matter physics 131-36. 39 -41] where 2D semiconducting thin 
films and ID semiconducting nanowires with proximity in- 
duced i-wave superconductivity and a Zeeman field have been 
shown to support Majorana fermions in order parameter de- 
fects such as vortices and boundaries. Such localized Majo- 
rana fermion excitations are thought to support non-Abelian 
statistics which has been proposed as a key resource for fault- 
tolerant topological quantum computation (TQC) li42[l . 2D ul- 
tra cold fermionic gases with laser induced Rashba SOC and 
a Zeeman field in the TSo phase can also support Majorana 
fermions in vortices and sample boundaries and thus can po- 
tentially act as a platform for fault tolerant TQC. In contrast, 
for 3D Fermi superfluid as in Fig. [T](b), the blue line deter- 
mined by the Chandrasekhar-Clogston like condition F = A 
is a topological phase boundary between fully gapped NTS 
(F < A) and gapless TS (F > A) phases. In addition, within 
the gapless TS phase, Fc is a quantum critical point of the 
topological phase transition between the TS phases with two 
(TS2) and four (TS4) nodes. The fermions with linear disper- 
sion relations near each node in the TS2 and TS4 phases are 
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a manifestation of the relativistic Weyl fermions. The bulk 
nodal excitation spectra of these phases are also topologically 
protected and a topological invariant characterizing the nodal 
points (in the context of He3) was formulated by Volovik |37]. 
Recently, similar topologically protected Weyl fermion modes 
have been shown to exist also in condensed matter systems 
such as ferromagnetic superconductors llssll . 



B. Chemical potential 

Chemical potential can be calculated by solving the order 
parameter equation with fixing the total number of the system 
self-consistently. With Zeeman field F and interaction v as 
tuning parameters, one can have chemical potential as a func- 
tion of r and V, that is, fi = fi(r, v). 

In Fig. |2l we show the chemical potential as a function of 
F for a given interaction v and SOC aKf. Fig. |2](a) and (b) 
are plots of the chemical potential fi and its derivative with re- 
spect to F , ^, for two dimensional Fermi gas at E/, = O.SEf 
and aKf = 0.4. In the side of low Zeeman field F, where 
F < Ff = ^J|J^ + and the system is in the non-topological 
superfluid phase (NTS), is slowly varying as F increases. In 
the region of strong Zeeman field, where F > F^. and the sys- 
tem is in a gapped TSo phase, fi decreases fast with increasing 
F. Near the quantum critical point F^., the slope changes dis- 
continuously as can be seen clearly in the plot of ^ in Fig.|2] 
(b), indicating the change in the curvature. 

As for the 3D Fermi superfluids, topological phase transi- 
tion from a fully gapped NTS to a gapless TS4 state occurs at 
Zeeman field F = Fa = A, and TS4 to TS2 phase transition 
occurs at F = F^. - -sjfi^ + A^. Fig.|2](c) and (d) illustrate fi as 
a function of F at unitarity (l/KfOs - 0) with the same SOC 
strength as that in 2D fermion gas. Even though the signa- 
tures of topological phase transitions across the phase bound- 
aries are weaker than the one for the 2D Fermi system, fi in 
the fully gapped NTS phase (F < Fa) increases with increas- 
ing F, while fi decreases with increasing F when it is in the 
TS2 phase (F > F^.). The curvature of the chemical potential 
versus the Zeeman field changes as the system moves from 
the gapped to the gapless phases and vice versa. In contrast 
to the 2D Fermi superfluid, however, ^ plot doesn't exhibit a 
cusp at the phase boundary because of the existence of the TS4 
phase surrounded by the gapped NTS and gapless TS2 phases. 
It is when ju = that the two distinct critical points Fa and F^ 
merge (Fa = Fc) and a tricritical point occurs, leading to the 
vanishing of the TS4 phase and the emergence of the cusp in 



^ even in a 3D Fermi superfluid system. 



C. Compressibility 

The compressibility can be obtained from the knowledge 
of the pressure, as defined by, kt - -il/V){dp/dV). But in 
the grand-canonical ensemble, where we need to fix the aver- 
age number of particles, the isothermal compressibility can be 
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FIG. 2. (Color online) Plot of chemical potential and its derivative 
with respect to F, ^ for (a)-(b) 2D Fermi gas at Ei, = 0.5Ef and 
(c)-(d) 3D Fermi gas at l/Kpa, = with aKp = 0.4. Blue and 
red crosses indicate the critical points Fa = A and F^ = -sjff + A-, 
respectively. One finds that fully gapped NTS phase occupies at F < 
Fc for 2D and F < Fa for 3D gases, respectively. 



directly written as. 



_}_(dN\ 



(11) 



In the superfluid state, we need to consider dependence of 
A. From the order parameter equation (5Q/5A) - 0, we have 
(d/dix)(dn/dA) = 0, leading to 



\df^lT,v 



ldN_ 



(12) 



Thus the isothermal compressibility in the superfluid phase 
can be written as 



Kt 



1 



dN\ 



(13) 



In Fig. [3] we present the isothermal compressibility kt as a 
function of F for a given interaction v at T = with aKf = 0.4 
to investigate the signature of the topological phase transi- 
tions. Since kj is related to the first derivative of the total 
number with respect to ji, that is, kt ~ '-^It.v, it is expected 
to show the signatures of the topological phase transitions as 
seen in the plots of ju in Fig.|2l In Fig. [3] kt shows a decreasing 
trend in the superfluid phase with increasing F, while it is an 
increasing function of F in normal phase as will be shown in 
Fig. |5] For 2D uniform Fermi superfluids, at the critical Zee- 
man field Fc, a sharp peak emerges as shown in Fig.|3](a). This 
peak in the compressibility corresponds to the phase boundary 
between the 2D gapped NTS and the 2D gapped TSo phases. 
In our calculations we find that the compressibility provides 
the strongest signature of the phase transition to the 2D topo- 
logical phase with the Majorana fermion excitations. In 2D at 
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FIG. 3. (Color online) Plot of isothermal compressibility kt as a 
function of T for (a) 2D Fermi gas at Ej, = Q.SEp and (b) 3D Fermi 
gas at 1/A'f a, = with aKp = 0.4. Blue and red crosses indicate the 
critical points Fa = A and = -y'p^ + A^, respectively, same as in 
Fig.IHand Fig.|2] 



zero temperature there is no additional phase transition from 
the TS() phase to the normal phase at higher values of F. This 
is because, with a finite Rashba SOC, the superfluid pair po- 
tential always remains non-zero even at large values of the 
Zeeman field. When the temperature is included as another 
parameter on the phase diagram, a normal phase with A = 
emerges even in 2D (see below). 

As for the 3D case. Fig. |3] (b) shows that KriT) again 
decreases with increasing F in the superfluid phase, and 
shows the signature of the topological phase transitions by a 
smoothly peaked curve near F = F^. As already seen in the 
plot of yU vs. F, Kt of 3D fermion system in the TS4 phase, 
which is sandwiched between the gapped NTS and gapless 
TS2 phases, is a fast varying function of F, while the com- 
pressibility in the gapped NTS phase is slowly varying. 



IV. THERMODYNAMIC SIGANTURES (T + 0) 

Experimentally, the thermodynamic properties of a given 
substance is determined by measuring an equation of state 
(EoS), such as the pressure piji, T) as a function of the chem- 
ical potential /i and the temperature T . Equivalently, replac- 
ing the pressure by the density n — (^)^, one can determine 
the density EoS, n(yU, T). The local gas density n{V) can be 
measured 126*] as a function of the local trapping potential V 
from in situ absorption images of a trapped strongly interact- 
ing Fermi gas at a Feshbach resonance. From the definition 
of the compressibility kt - ^ chemical potential 

as a function of temperature T is attainable experimentally 
without numerical derivatives of the data ll26ll . Then, from 
the energy, pressure, and the chemical potential, the entropy 
S - j{E + pV - fiN) and the specific heat cy can be obtained 
as a function of temperature T. 



A. Thermodynamic Properties 

From the mean-field grand partition function 

Z = Tr(e-(«-''^>/^) (14) 



we can calculate the corresponding grand potential function 
Q = -rinZ = (H) -TS - p{N) = -pV, where S is the 
entropy and p is the pressure of the system. The differential 
of the grand potential function is 



dQ. ^-SdT- pdV - {N}dfj., 



(15) 



allowing us to obtain interesting thermodynamic properties. 
The explicit expression of Q. in terms of the Hamiltonian un- 
der consideration is 



Then, the pressure of the ultracold Fermi gas can be ob- 
tained through the grand potential function via p = -Q./V. In 
the thermodynamic limit, it reads 

n = > -ln{l+e '^■■'')+ I ^r^k + — 

^ 2^ J (2nr ^ ' J (271)0^" g 

(17) 



The entropy can be obtained by 



(18) 



Then, the heat capacity can be obtained from the knowledge 
of the entropy. 



Cv 



dS_ 

dT 



(19) 



In the superfluid phase, A depends on temperature T as 



dS\ ld^Q.V 



(20) 



(21) 



leading to the specific heat in the superfluid phase 
-i—t (—]' 

B. Finite Temperature Pliase Diagram 



Fig. |4]presents mean-field (T-T) phase diagram for 2D and 
3D homogeneous Fermi gas in the presence of Rashba-type of 
spin-orbit coupling. We choose the interaction parameter v as 
a two-body binding energy Ei, - 0.5Ef for 2D Fermi gas (a), 
while 1 /KfO, = for 3D Fermi gas (b) with Rashba SOC 
aKf - 0.4. With increasing temperature T, the superfluid 
order parameter A decreases and undergoes the second order 
phase transition to normal phase (A = 0) denoted as N at T = 
Tc- 

In Fig. m (a), within the superfluid phase of 2D fermion 
gas, topological phase transition between the gapped NTS 
phase (F < Fc) and the gapped TSo phase (F > Fc) occurs 
atT - Ti <Tc indicated with a red line. In addition, the TSo 



6 



(a) (b) 
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FIG. 4. (Color online) Plot of (F - T) phase diagrams for (a) 2D 
Fermi gas at Et = O.SEf and (b) 3D Fermi gas at 1 IKfO^ = with 
aKf = 0.4. Red and blue lines denote F^. = ^jJjpT^ and Fa = A, 
respectively. Black lines represent Fq above which the order param- 
eter vanishes, or A = 0. The tricritical temperatures, at which F^ 
meets with Fo (black circles), for 2D and 3D gases are T - 0.43£/r at 
F,. = Fo = Q.(>2Ef and T = QMEp atF,. = Fo = 0.42£f, respectively. 



phase vanishes at the tricritical temperature - - Q.A3Ef 
with our choice of parameters Eb and a. In Fig.|4](b), T2 (blue 
line) determined by the Chandrasekhar-Clogston-like condi- 
tion, r = A, is a phase boundary between the NTS phase 
(T < T2) and the TS4 phase (T2 < T < min[rc,ri]). It is 
when r = that T2 - Tc occurs and the system above 
becomes the normal state (A - 0). Thus, for F = 0, there is 
only one finite temperature phase transition at Tc{= T2) which 
is a transition between the non-topological superfluid (NTS) 
(T < Tc) and the normal state (T > Tc). In a similar man- 
ner as in the 2D system, TS4 phase vanishes at finite F, where 
Tc = Ti= 0A3Ef at unitarity (1 /Kpa, = 0). 



C. Thermodynamic Quantities as a function of temperature 

Since we learned that the chemical potential and the isother- 
mal compressibility at zero temperature can show the evi- 
dence of the topological phase transitions as a function of the 
tuning parameters such as the Zeeman field F and the interac- 
tion parameter v, we first consider below the chemical poten- 
tial and the compressibility at finite temperatures. 

In Fig. |5] we show the chemical potential and isothermal 
compressibility as a function of the reduced Zeeman field 
F/Fo for 2D on the left column and 3D on the right column 
at both T > T[ (dashed) and T < Ti (solid). We choose 
the Rashba SOC aKf - 0.4 and the interaction parameters 
El, - Q.SEf for 2D and l/Kfa, = for 3D Fermi gases, 
respectively. Fig. |5] (a) and (b) show that, at low tempera- 
tures, T < T\, kt for both 2D and 3D fermions exhibit peak 
at or near the topological critical point, and a discontinuity 
at the phase transition to the normal state at F - Fq. Mean- 
while, at high temperatures (dashed), T > Ti , it is hard to see 
the signature of the topological phase transition between the 
gapped NTS and gapless TS4 phases. We notice that as tem- 
perature increases kj increases in the superfluid phase, while 
it decreases in the normal phase. Fig. |5](c) and (d) are plots 
of both 2D and 3D chemical potentials as a function of F/Fq 




r/Fo r/To 

FIG. 5. (Color online) Plot of isothermal compressibility kj and 
chemical potential as a function of F/F^. for (a),(c) 2D and (b),(d) 
3D Fermi gas at temperature T = Q.lEf (solid) and T = 0.45£'f 
(dashed). We indicates topological phase transitions with red (TS4 
to TS2) and blue (NTS to TS4) crosses within the uniform superfluid 
phases. One finds that normal (N) to superfluid phase transitions oc- 
cur at Fo and are characterized by the discontinuity of the isothermal 
compressibility kt, while the chemical potential fi remains a contin- 
uous function of F as well as a function of T. 



for T - O.lEp (solid) and T — OASEp (dashed), respec- 
tively. It shows that fi is an increasing function of T and /u 
in the superfluid phase is varying slowly compared with that 
in normal phase. Even though the signatures of the topologi- 
cal transition, which are strong at zero temperature, are we ak- 
ened by the thermal fluctuations, the characteristic features for 
the variations of the thermodynamic quantities near the phase 
transitions remain the same. 

In Fig. |6l we present the other relevant thermodynamic 
quantities in 2D and 3D Fermi gas, such as pressure p, entropy 
S , and specific heat cy, which can be observed expeiimentally, 
as a function of F/Fo at T = O.lEf (solid) and T - QA5Ef 
(dashed). We choose Eb = Q.SEf for 2D and l/Kpas - for 
3D Fermi system with Rashba SOC aKf = 0.4. The main 
conclusion from these calculations is that the signatures of the 
topological phase transitions are much weaker in these quan- 
tities than those in the case of ordinary phase transitions from 
the ordinary i-wave superfluid to normal phases. 



V. CONCLUSIONS 

In summary, we have investigated the two- and three- 
dimensional uniform fermion ^-wave superfluids in the pres- 
ence of Rashba spin-orbit coupling together with Zeeman field 
as they undergo topological phase transitions as a function of 
the interaction parameter v and Zeeman field F at zero tem- 
perature and as a function of Zeeman field F at finite temper- 
atures. In 2D 5-wave Fermi gas, for F = Fc = -y/Ju^'+'A^ the 
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FIG. 6. (Color online) Plot of pressure p, entropy S , and specific 
heat Cv as a function of T/Tq, where Fq corresponds to the Zeeman 
field at A = 0, for (a),(c),(d) 2D Fermi gas at £(, = 0.5^^ and 
(b),(d),(f) 3D Fei-mi gas at l/Kfa, = at T = O.IEf (solid) and 
T = OASEf (dashed). We used the red and blue crosses for Fj. and 
Fa, respectively. One notices that the signatures of the normal (N) 
to superfluid phase transitions are visible in both 2D and 3D Fermi 
gases, characterized by the suppression of p and 5 in the superfluid 
state and the discontinuity of cy. On the other hand, it is difficult to 
observe the topological phase transitions even with some weak sig- 
natures. 



superfluid phase is identified wifli a gapless topological super- 
fluid (with one node at A: = 0), while it becomes a fully gapped 
non-topological superfluid phase at F < Fc and a fully gapped 
topological superfluid phase TSo, capable of supporting Ma- 
jorana fermion excitations at order parameter defects such as 
vortices, at F > Fc. The TSo phase has recently received 
wide attention in condensed matter physics ||3T1 - |36[ [39I - I4TI1 
where 2D semiconducting thin films and ID semiconduct- 
ing nanowires with proximity induced i-wave superconduc- 
tivity and a Zeeman field have been shown to support Ma- 



jorana fermions in order parameter defects such as vortices 
and boundaries. Such localized Majorana fermion excitations 
are thought to support non-Abelian statistics which has been 
proposed as a key resource for fault-tolerant topological quan- 
tum computation (TQC) [42.1 . In contrast, 3D i-wave Fermi 
gas is a fully gapped non-topological superfluid at F < A, 
becomes a gapless topological superfluid with four topolog- 
ically protected Weyl fermion points iIstIi at A < F < F^ 
and another topological superfluid with two topologically pro- 
tected Weyl fermion points at F > F^. With our calculated 
mean-field phase diagrams in both 2D and 3D we investigated 
the thermodynamic properties across the relevant topologi- 
cal phase transitions, including the isothermal compressibil- 
ity, the chemical potential, the pressure, the entropy, and the 
specific heat as a function of the Zeeman field F and the inter- 
action parameter v. In particular, the isothermal compressibil- 
ity of the system supporting the Majorana fermions, which is 
the 2D fermion gas with Rashba SOC together and a perpen- 
dicular Zeeman field, reveals the topological phase transition 
via a strong peak at the critical point F - F^ (Fig. O. Fur- 
thermore, the phase boundary can also be identified by the 
chemical potential through the change of the curvature as a 
function of F (Fig.|2|. For the 3D Fermi gas, the topological 
phase transition from the fully gapped non-topological phase 
to the gapless topological phases can be probed by the posi- 
tion of the peak in the isothermal compressibility as a func- 
tion of F, even if the signatures of the topological transitions 
are weaker than those in the 2D system. The signatures of the 
topological phase transitions at low temperature disappear as 
the temperature goes above the critical temperature Ti. The 
other thermodynamic quantities, such as the pressure, the en- 
tropy, and the specific heat, have only weak signatures of the 
topological phase transitions (that is, the transition between 
a conventional (non-topological) superfluid and a topological 
superfluid such as the Majorana and the Weyl phases), while 
these are useful quantities to identify the phase transitions 
from the normal phase to the topological superfluid phases 
with decreasing values of the Zeeman field. 
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